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Abstract 

We consider Af = 4 supersymmetric Yang-Mills theory with SU(N) gauge group 
at large N and at finite temperature on a spatial S 3 . We show that, at finite weak 't 
Hooft coupling, the theory is naturally described as a two dimensional Coulomb gas of 
complex eigenvalues of the Polyakov-Maldacena loop, valued on the cylinder. In the 
low temperature confined phase the eigenvalues condense onto a strip encircling the 
cylinder, while the high temperature deconfined phase is characterised by an ellipsoidal 
droplet of eigenvalues. 



1 Introduction 



1.1 Known results on the phase structure of thermal Af = 4 SYM 

At present N = 4 supersymmetric Yang-Mills (SYM) theory provides the canonical ex- 
ample of an explicit gauge theory/string theory duality at large N [1]. The extension of 
the correspondence to finite temperature [2] suggests that the duality may allow the phase 
structure of the theory to be fully mapped out as a function of temperature and coupling. 
From a gravitational perspective, the finite temperature theory furnishes a unique window 
into black hole physics. A key question, for instance, is whether stringy corrections funda- 
mentally alter the nature of black holes in strongly curved spacetimes, or whether weakly 
coupled plasmas are qualitatively similar to black holes. 

The determination of the phase structure of N = 4 SYM theory at finite temperature 
on a spatial S* 3 began with Witten's identification [2, 3] of the first order Hawking-Page 
transition in the bulk as a strong coupling deconfinement transition. It was later noticed 
that a similar first order transition [4, 5] occurs at a critical temperature in the free field 
theory on S* 3 at large N. In both of these transitions, at zero and strong coupling, the 
free energy of the system jumps from 0(1) to 0{N 2 ) at the critical temperature, as the 
theory deconfines. This suggests that the phase transitions should be identified via a phase 
boundary running from zero to strong coupling. In fact this boundary turned out to have 
some structure of its own, as is described in [5]. 

The deconfinement transition at weak coupling was given an elegant treatment in [5], 
where the theory was systematically shown to reduce to a unitary matrix model for the 
Polyakov loop. This is the holonomy of the SU{N) gauge field around the Euclidean time 
circle of the thermal field theory 

U = Pe i f AodT . (1) 

In the low temperature confining phase (TvU k ) = 0, while at high temperatures these 
traces are nonvanishing. Such behaviour of these traces implies that, in the large N limit, 
the distribution of eigenvalues of the Polyakov loop is uniform in the confining phase and 
nonuniform at high temperatures. 

1.2 The role of scalar fields 

This work was motivated by two open questions from previous treatments: 

• The order parameter for the weak coupling deconfinement transition is the eigenvalue 
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distribution of the Polyakov loop. However, the order parameter which can be com- 
puted at strong coupling, and which is natural from the viewpoint of the dual string 
theory, is the Polyakov-Maldacena loop [6]. To compare the phases at strong and 
weak coupling, we would like to use the same order parameter. 

• At zero 't Hooft coupling all fields except the zero mode of Aq (or equivalently, the 
Polyakov loop), are massive and can be integrated out. This is possible because all 
the scalars have positive conformal mass term. However, we shall see that quantum 
effects at nonzero 't Hooft coupling A cause condensation of the scalar fields as well 
as A) [7]. 

These two questions are related. Recall that the Polyakov-Maldacena loop [6] is 

W = p e if[A o +i<S>j0 J (r)]dT ^ £ 2 ) 

where 6 J is a unit vector in M 6 . We can see that if Aq and <&j have commuting nonvanishing 
vacuum expectation values, then the eigenvalue distribution of iAq + $j is well defined and 
will be directly related to the computation of (tvW k ). 

In this letter we study the eigenvalue distribution of <£j + iAq at finite weak coupling 
as a function of temperature. This is a distribution on the complex cylinder IxS 1 . The 
resulting two dimensional framework is the necessary generalisation to nonzero coupling 
of studies of the free theory, such as [5], which were phrased in terms of the eigenvalue 
distribution of Aq on S 1 . 

In the limits of low and high temperature the eigenvalue distribution is given precisely by 
a 2D Coulomb gas in an external potential. At intermediate temperatures, the interaction 
potential has deviations from the Coulomb law. We compute the distributions analytically 
in these limits: the results are constant density droplets that are shown in figures 1 and 2 
below. 

2 The two dimensional Coulomb gas 
2.1 One loop effective potential 

To understand the phase structure of the theory at finite, weak 't Hooft coupling we need 
the quantum effective potential generated at finite temperature. Such an effective potential 
can be a function of all homogeneous modes in the theory on S 3 x S 1 . In particular it will 
depend on the zero modes of the six scalar fields of the N = 4 theory and the Polyakov loop 
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(or Aq). Fortunately, this potential has very recently been computed at one loop order [7]. 
Several consistency checks were also performed in that paper. There are three contributions: 
the classical term and then the bosonic and fermionic fluctuation determinants obtained by 
integrating out the respective Kaluza-Klein harmonics on x S 1 . The potential turns out 
to be a gauge invariant function depending only on the eigenvalues {4>Jp} of the adjoint 
scalars and those of the time component of the gauge field {0 P } 

SeB[<p p , P ] = S^ + S^+S^. (3) 

The classical term is 

S^=f3R^ (4) 

p=0 

which is just the scalar mass term due to conformal coupling to the curvature of the S* 3 
of radius R. Furthermore, A = gyu-^ "C 1 is the 't Hooft coupling and T = jjj, the 
temperature. The reason it is possible to have simultaneously diagonal scalar and A 
background values can be traced to the Af = 4 classical potential and the fact that it allows 
a Coulomb branch moduli space on R 3 ' 1 . In the presence of diagonal background fields, the 
off-diagonal fluctuations are massive and appear quadratically in the action. This allows 
them to be consistently integrated out with vanishing expectation values. 

Integrating out the gauge fields, scalar and ghost fluctuations yields the bosonic contri- 
bution to the one loop potential 

JV-1 

p,q=0 

oo 

+ ^(2^ + 3)(2^ + l)log 
e=o 

and 

JV-1 oo 

<? f (1) = {pct( Vpq )-J2u(e + i)iog 

p,q=0 l=\ 

In these expressions pq = 6 P — 6 q and 

j=i j=i 

The functions and Cf are the bosonic and fermionic Casimir energy contributions in 
the presence of background expectation values for the adjoint scalar fields [7]. Their sum 



cosh/3-, /ip 2 



cos 0. 



pq 



.log 2 



+ c.c 



(5) 



1 + e 



-^(t+\) 2 R- 2 +V 2 Pq +iO P 



+ C.C 



(6) 
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yields the zero temperature one loop effective potential of N = 4 theory on 5 3 . For present 
purposes it will suffice to know that this potential has the expansion 



Cb{<£pq) + Cf{ip pq ) = — 



A + (^)2(io g 2-i)+0(i?V^) + --- 



(8) 



near the origin. 

The problem now is to find the eigenvalue distribution that minimises the effective 
potential (3). 

2.2 Instability of the ring configuration and short distance repulsion 

Previous works have considered condensation of the gauge field eigenvalues 6 P while assum- 
ing that the scalar eigenvalues remained zero: 4>j p = 0. We will instead look for general 
distributions in which both sets of eigenvalues can condense. Let us introduce the complex 
field 

Zp = ^±^L. (9) 
The z p live on a cylinder: 9ie z p ranges from — oo to +oo but 3m z p is periodic with range 

7T. 

In our definition of the complex variables z v and in what follows we will omit 50(6) 
vector indices on cf> j p and think of it as the radial mode in 50(6) space. A non-zero radial 
mode implies symmetry breaking. This is possible on 5 3 x 5 1 due to the large N limit which 
plays the role of the thermodynamic limit [7]. Hence we set ip p = <p p and ip pq = 4> p — <p q . 

The first statement we can make is that although (f) p = is always a solution, call it the 
ring configuration, it is never a stable configuration at finite 't Hooft coupling. This is due 
to a strong repulsive force between the eigenvalues at short distances, originating from the 
first line of (5). The term may be rewritten as 

i (log 2 + log cosh PyJ^f q - cos pq ) = log 2 + log I sinh z pq \ . (10) 

This term reduces to the familiar log sin when (j) p = 0. As z p ^ z q we see that the eigen- 
values experience a 2D Coulombic repulsion with potential — log \ z p — z q \. At large N there 
is no other term in the action (3) that can compete with this term at short distances and 
therefore the eigenvalues always spread out into two dimensions. Let us see this explicitly. 
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2.3 Low temperature limit 



Consider first the low temperature limit of the confining phase. That is, take ^ <C 1. The 
action (3) at low temperatures becomes 



Srt<ci 



3/3 
16R 



log 2 



+ 



Nn 2 R 

(3\ 



p=0 



N-l 

log | sinh z, 

p,q=0 



/«/ 



(11) 



As well as the terms in the action which are exponentially suppressed at low temperatures, 
we have also dropped all the Casimir energy terms from (5) and (6). This will be valid if 
R<t>Jp ^ 1- We will see a posteriori that this is true for the solution we are about to present. 

The action in (11) is an elegant reformulation of the problem as a Coulomb gas in two 
dimensions with an external potential. This system is closely related to normal matrix 
models. The interaction potential log | sinh z pq | is precisely the Coulomb potential on a 
cylinder. 

In writing down the equations of motion following from (11) it is convenient to pass to 
the continuum large N limit. This is done by introducing an eigenvalue distribution 

N-i , 

J2^ N d 2 zp(z,z). (12) 



q=0 



The continuum limit implies the normalisation condition 1 

d 2 zp(z, z) = 1 , 



while the equations of motion become 

(3\ 



z + z 



2it 2 R 



d 2 z'p{z' , z') coth(z 



(13) 



(14) 



As is standard for these types of system, we can immediately derive from (14) that the 
eigenvalue distribution is everywhere either a fixed constant, or it is zero. This follows from 
acting on (14) with 8, and gives 



p(z,z) 



f if zeB 



(15) 



if z i B 

Here B is a subset of the complex cylinder. In other words, the eigenvalue distribution 
takes the form of a uniform density 'droplet' on the cylinder. The equations of motion (14) 
therefore imply that we need to find a domain B such that 



+ z = — d 2 z' coth(z - z') . 
27T J B 



(16) 



1- We work with the convention that d 2 z = 2dxdy, where z = x + iy. 
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From symmetry considerations it is not difficult to see what form the droplet has to 
take. The external force acting on the eigenvalues acts parallel to the real axis and pushes 
the eigenvalues towards the imaginary axis. The repulsive force between the eigenvalues 
must balance this effect to lead to a finite width band. This is illustrated in figure 1 below. 
We will show momentarily that a band of constant width 2 A does indeed solve (16), for all 
A. Given this fact, the normalisation condition (13) determines the width A of the band to 
be 



n 




Figure 1 : The low temperature distribution of eigenvalues in a band of width 2 A = 2 ^r 
on the z cylinder. 

The equation of motion (16) is checked by doing the double integral. Set z' = x + iy so 
that the integral becomes 

f-A 



— I d 2 z' coth(z — z') = — f dx j dycoih(z — x — iy) 
27r Jb K J -A Jo 

f(z+z)/2 f 

= dx- 

J-A J(; 



(z+z)/2 r A 

dx 

(z+z)/2 

= z + z. (18) 



In going from the first to the second line here, one needs to be careful to keep track of the 
branch cut in logsinh. 

Given the band solution, we can justify the low temperature approximations made to 
obtain the action (11). We see that R<p ~A<1. It follows that evaluated on the solution, 
the Casimir energy terms are suppressed by a power of the 't Hooft coupling relative to 
the terms kept in (11). It furthermore follows that because the spread RA(f> ~ A makes 
no reference to the temperature (3, the condensate has a finite zero temperature limit. 
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Consistent with confinement, and with strong coupling results, the band solution yields 
(TvW k ) = (TrU k ) = due to cylindrical symmetry of the distribution which causes all 
phases to average out to zero. 

Having established that the finite band solves the equations of motion, we can check 
that indeed the configuration has lower action than the ring z + z = 0, which we argued 
was unstable. Evaluating the action (11) on the solution in the continuum limit gives 

3 A 



'-'band ^ 



16 12vr 2 



(19) 



The ring solution at low temperatures is uniformly distributed along the ring, 9 P = 2-Kp/N, 
which gives the action 

>Sring = — > iSband • (20) 

2.4 Deconfinement transition and high temperature limit 

If we raise the temperature and consider the effects of the exponentially suppressed terms 
in (5) and (6), then we find that the constant density band solution we have just described 
remains a solution. This occurs because the action may be expanded in terms of cos(n0 pq ) 
factors, which give zero when summed over the uniform distribution on the 9 circle. 

The behaviour of these sums of exponentials is identical to that described in [5] for the 
Polyakov loop matrix model. Let us recall from that work that when (3/R becomes of order 
one, then although the uniform distribution in the 9 direction remains a solution, it is no 
longer the absolute minimum. This is the deconfinement temperature. The presence of 
R(fpq in the exponent, the new ingredient in our computations, will not change this story 
qualitatively. It will however play a role in determining the order of the transition, c.f. [5]. 

Above the deconfinement transition, the force between eigenvalues in the circle direction 
becomes attractive and the eigenvalue distribution develops a gap along the 9 circle. We 
postpone a careful discussion of the behaviour of eigenvalue droplets near the transition for 
future work [9]. The story is complicated by likely non uniformities in the 4> direction. 

A regime that is again analytically tractable is the high temperature limit, ^ S> 1. From 
previous work [5] , our intuition is that at high temperatures the eigenvalues will clump into 
a small droplet. Let us take, simultaneously to the high temperature limit, the limit in 
which 9 p , (5(pq <C 1 in the action (3). Similarly to before, we will verify a posteriori that this 
condition is satisfied by the solution. It is convenient to again introduce the complex z p as 
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we did in (9). The action in these limits becomes 

Ntt 2 R N ~ 1 N ~ 1 ir 2 R 3 N ~ 1 

Str»i = — ^— Y, ( z p + h? ~ E lo S M + ~7p~ Yl ~ z pi + 6 I^H • ( 21 ) 

P p=0 p,q=0 P p,q=0 

Although the first two terms in this expression are naively subleading, they become impor- 
tant depending on the 't Hooft coupling and proximity of the eigenvalues. We now take the 
continuum limit as we did at low temperatures (12). The same arguments as before imply 
that the eigenvalues form a droplet B of constant density 

6R 2 Xl 



ttR 



1 + 



P 2 



if z G B. (22) 



The equation of motion becomes 



1 



tanhfiz + z = / d 2 z' — ^— , (23) 

where we defined 



Itt J d z — z' ' 



1 - 



B 
2R 2 X 



ta ^ = T-4' ( 24 ) 

1 + ~W 

The equation (23) has been simplified by the observation that if the symmetry z — > — z is 
unbroken, then we must have z q = J2 q % = on the solution. 

Before solving the equation of motion, the normalisation condition (13) tells us that the 
area of the droplet on the complex plane will be 

Area < B > = i^T7V- < 25 > 

This area goes to zero as (3/R — > at high temperatures, and therefore our initial assumption 
on (j) p and 6 P was justified. 

The droplet shape for the Coulomb gas problem (23) is known to be an ellipse [8]. This 
is illustrated in figure 2 below. The semi-axes of the ellipse are of length re _M along the 
real axis and re^ along the imaginary axis, where [i is given in (24) and irr 2 equals the 
area (25). The eccentricity is therefore controlled by the ratio XR 2 /P 2 . Interestingly, it is 
precisely this ratio which controls the validity of perturbation theory at finite temperature 
on S 3 . The one loop effective potential and results we have deduced from it are valid for 
all temperatures T < [7]. 



3 Discussion 

In this letter we have shown that weakly coupled N = 4 SYM theory at finite temperature on 
a spatial S 3 is described at large TV" by an eigenvalue distribution on the complex cylinder. 
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Figure 2: The high temperature distribution of eigenvalues in an ellipse (not to scale). 

Remarkably, these eigenvalues are naturally associated to the Polyakov-Maldacena loop, 
the order parameter known to emerge at strong coupling. At high and low temperatures 
we determined the distribution by solving an associated 2D Coulomb gas problem. The 
deconfinement transition of the theory on S 3 corresponds to a 'band-to-droplet' transition 
on the cylinder. A key point we emphasised is that the eigenvalues of the scalars as well as 
Aq condense. Outstanding questions are the implications of scalar condensates for the dual 
string theory and the relevance of the 'band-to-droplet' picture for the emergent Big Black 
Hole phase at high temperatures. 

It is of immediate interest to understand the droplet picture in the vicinity of the phase 
transition, and also whether the strongly coupled theory may be formulated in the Coulomb 
gas language. This could lead to a new approach for counting black hole microstates. Some 
progress studying higher traces of the Polyakov-Maldacena loops at strong coupling was 
recently made [10]. The Coulomb gas framework we have presented has a description in 
terms of fermions, suggesting a connection to a different droplet picture of certain BPS 
states recently developed for the zero temperature theory [11, 12]. 
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